o 



(N 



Photometric transit search for planets around cool stars from the 
Western Italian Alps: a site characterization study ^ 

M. Damasso^'^, P. Giacobbe^, P. Calcidese^, A. Sozzetti^, M.G. Lattanzi^, A. Bernagozzi^, 

E. Bertolini^, and R.L. Smart^ 



ABSTRACT 

We present the results of a site characterization study carried out at the Astro- 
nomical Observatory of the Autonomous Region of the Aosta Valley (OAVdA), in 

the Western Italian Alps, aimed at establishing its potential to host a photomet- 

Cm ' 

rjj ■ ric transit search for small-size planets around a statistically significant sample of 

(-H ' nearby cool M dwarfs. For the purpose of the site testing campaign, we gathered 

Qh' photometric and seeing measurements utilizing different instruments available at 

O ■ the site. As in any search for new locations for astronomical observations, we 

'^ • gauged site-dependent observing conditions such as night-sky brightness, pho- 

^ ■ tometric precision, and seeing properties. Public meteorological data were also 

used in order to help in the determination of the actual number of useful ob- 
^ ' serving nights per year. The measured zenithal \^-band night-sky brightness is 

typical of that of very good, very dark observing sites. The extinction registered 



CN ■ at V band is not dissimilar from that of other sites. The median seeing over the 

period of in situ observations is found to be ~ 1.7". Given the limited duration 

(3 ■ of the observations, we did not probe any possible seeing seasonal patterns, or 

5 ■ the details of its possible dependence on other meteorological parameters, such as 

wind speed and direction. Moreover, our data show that the seeing at the obser- 

i^ ■ vatory was reasonably stable during most of the nights. The fraction of fully clear 

'}_j ■ nights per year amounts to 39%, while the total of useful nights increases to 57% 

assuming a (conservative) cloud cover of not more than 50% of the night. Based 

on the analysis of photometric data collected over the period May- August 2009 

for three stellar fields centered on the transiting planet hosts WASP-3, HAT-P-7, 
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and Gliese 436, we achieve seeing-independent best-case photometric precision 
(jph < 3 nimag (rms) in several nights for bright stars {R < 11 mag). A median 
performance Uph ~ 6 mmag during the observing period is obtained for stars 
with i? < 13 mag. A byproduct of the significant amount of photometric data 
collected in the stellar fields of WASP-3 and HAT-P-7 was the identification of a 
handful of new variable stars, four of which were presented and discussed here. 
Our results demonstrate the OAVdA site is well-poised to conduct an upcoming 
long-term photometric survey for transiting low-mass, small-size planets around 
a well-defined sample of M dwarfs in the solar neighborhood. 

Subject headings: Site testing - stars: individual: WASP-3, GJ 436, HAT-P- 
7, GSC2.3 N208000215, GSC2.3 N2JH035417, GSC2.3 N2JH000428, GSC2.3 
N2JH066192 - planetary systems - Stars: variables: general - Techniques: pho- 
tometric 



1. Introduction 

The discovery by Charbonneau et al. (2009) of a ~ 6.5 M® planet (a 'Super Earth') 
with a radius of ~ 2.7 R^ in transit across the disk of the late-M dwarf GJ 1214 with 
a period of ~ 1.6 days has heralded the first success for a novel approach to ground-based 
photometric searches for transiting planets. Until recently, two main typologies of search pro- 
grams for planetary transits had been adopted: 1) deep, narrow-field surveys using medium- 
to large-size telescopes with a narrow field of view (FOV) to monitor many faint sources, 
and 2) shallow, wide-field surveys utilizing small telescopes with a wide FOV to probe for 
photometric variability a large number of bright, relatively nearby solar- type stars. The 
former approach has produced so far the overwhelming majority of transiting planet detec- 
tions, and is the preferred choice for ongoing (CoRoT, Kepler) and planned (PLATO, TESS) 
space-borne transit searches (for a review, see Charbonneau et al. 2007). 

Among the many paradigms which the fast-developing field of extrasolar planet science is 
continuously reformulating, a special focus has been given in recent times to the choice of the 
stellar type to be probed for planets. Driven by the available instrument performance for a 
given technique (and by a small amount of hard-to-die Solar-System-centrism), initial planet 
search surveys (be they spectroscopic or photometric) focused exclusively on main- sequence. 



^ Based on observations made at the Astronomical Observatory of the Autonomous Region of the Aosta 
VaUey (Italy) 



-3- 



solar-type stars. This is changing, however, with increasing awareness of the need to extend 
the number of detected planets around varied star types, for a proper, global understanding of 
the complex processes of planet formation and evolution. With the primary goal of reaching 
sensitivity to short-period planets with (minimum) masses lower than the mass of Neptune, 
a number of dedicated radial- velocity surveys have been targeting low-mass stars (e.g., Endl 
et al. 2006; Zechmeister et al. 2009 and references therein). High-precision Doppler data 
(Rivera et al. 2005; Udry et al. 2007; Forveille et al. 2008; Mayor et al. 2009) have unveiled 
the existence of a population of low-mass, close-in planets with Mp sin i K.2 — 1Q M^ around 
cool starqj. At wider separations, Super-Earths around M dwarfs have also been revealed by 
microlensing surveys (Beaulieu et al. 2006). 

As for the transit technique, the combination of photometric data and Doppler mea- 
surements for targets at the bottom of the main sequence implies the possibility of deriving 
estimates of the bulk composition for Super-Earth planets (with typical radii in the range 
2 — 4 -R©), with requirements on the photometric precision necessary for their detection 
which are readily fulfilled from the ground (for details see for example Nutzman & Char- 
bonneau 2008). However, the specific choice of targets has imposed a re-definition of the 
observing strategies. Given the sparseness of bright M dwarfs throughout the whole sky, 
the wide-field approach can only be carried out effectively utilizing 4-m class telescopes, and 
shifting to near-infrared observing wavelengths in order to be able to probe tens of thou- 
sands of M dwarfs which would be otherwise too faint in the visible, as it is done for example 



by the WFCAM/UKIRT Transit Survey (WTS, see |http : //star .herts . ac .uk/RoPACSTI ). 
Alternatively, the narrow-field approach has been reformulated by designing transit search 
programs which aim at probing individually all the sparsely spread bright, nearby M dwarfs 
in the sky, utilizing a cluster of identical telescopes at one location, such as the MEarth 
project (Nutzman & Charbonneau 2008). The discovery by Charbonneau et al. (2009) is 
the proof of the cunningness of the latter method, given that only for transiting Super-Earth 
planets orbiting nearby small stars it is possible to carry out more in-depth characterization 
studies, such as spectroscopic investigations of their atmospheres. 

The ground-breaking discovery of GJ 1214b has, however, also highlighted some of 
the present limitations in our understanding of transiting systems composed of low-mass 
primaries and small-size planets. The first source of uncertainty comes from the existing 
discrepancies between theory and observations in the determination of the sizes of (late) 
M dwarfs which, as reported by Charbonneau et al. (2009, and references therein), are 
of the order of 10% — 15%. This problem, likely stemming from the need to introduce in 



^Evidence for their presence around solar-type stars has been presented by Bouchy et ah 2008; Fischer 
et ah 2008, and more recently by Howard et al. 2009; Vogt et al. 2009; Rivera et al. 2009 
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stellar evolutionary models the detailed treatment of the effects of non-zero magnetic fields 
on the properties of low-mass stars (e.g., MuUan & McDonald 2001; Ribas 2006; Torres et 
al. 2010, and references therein), has significant repercussions on the inferred composition 
of the planet. An important by-product of the long-term photometric monitoring of a large 
sample of low-mass M dwarfs is the discovery of new eclipsing binary systems with low- 
mass components (e.g., Irwin et al. 2009), leading to a sample of accurate measurements of 
their radii which can in turn provide a useful testing ground for improving stellar evolution 
models. Secondly, degeneracies in the models of the physical structures of Super-Earths, such 
as objects with very different compositions having similar masses and radii, also prevent one 
from exactly inferring their interior composition (Seager et al. 2007; Rogers & Seager 2010. 
See also discussion in Charbonneau et al. 2009) when only a mass and radius measurement 
are available. It is thus clear that more discoveries of transiting low-mass, small-radius 
planets over a range of periods around nearby, relatively bright late-type stars are highly 
desirable. This is of particular relevance as the further spectroscopic characterization studies 
of the atmospheres possible for these planets are one of the most effective ways to mitigate 
the degeneracies of exoplanet interior composition models (see for example Miller-Ricci et 
al. 2009). 

If the wild diversity of transiting giant planets properties is a proxy for what to expect 
from transiting Neptunes and Super Earths, and if the prediction from recent population 
synthesis models of planetary systems formation are on target, which envision increasing 
planet frequencies with decreasing planet mass, one then realizes how well justified is the 
effort being invested in present and planned transit search programs, both from the ground 
and in space. One lesson learned from the experience of wide-field ground-based transit 
surveys for giant planets around bright solar-type stars is that, in order to optimize the phase 
coverage and thus maximize the chances of success, a distributed network of telescopes should 
be built. Multi-site campaigns carried out by the HAT- Net and Super- WASP networks have 
contributed the large majority of the transiting planet discoveries to-date (the HAT-South 
network is now coming online, see Bakos et al. 2009). With this in mind, our group has 
started investigating the prospects for a Europe-based photometric transit search for small- 
size planets orbiting low-mass stars using a network of small telescopes. As a first step 
towards this goal, we have attempted to gauge the potential of the recently inaugurated 
(2003) Astronomical Observatory of the Autonomous Region of the Aosta Valley (OAVdA, 
[http : // www . oavda . it ) to become the primary node of such a network. Due to its short 
history, OAVdA is significantly involved in teaching and public outreach activities, but it 
also has an early record of collaborative efforts in a variety of scientific programs (e.g., 
monitoring of near-Earth objects and blazar events). The primary reasons for an in-depth 
characterization of the potential of the site for the purpose of a transit search program 
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are as follows. First, the Observatory is located in the North- West of the Italian Alps 
(45.7895 North, 7.47833 East), at 1,675 m above sea level, thus its high-altitude location 
and remoteness provide almost light-pollution-free night-sky conditions. Second, clear-sky 
observing conditions are realized at OAVdA on ~57% of the total nights in a year, which 
is a conservative estimate based on 4-yr statistics collected at the site (see § 14. 2p . Third, 
the OAVdA site has two platforms with electronically-controlled track-mounted ceilings that 
can be completely opened to the sky. One of these can host a variable number of telescopes, 
depending on system specifications and usable area on the platform. The availability of 
a single movable enclosure (the site is essentially unique in Europe) and of much of the 
associated support infrastructure at OAVdA will allow in perspective very substantial cost 
savings for the project. 

The generic design requirements for a ground-based photometric transit search for small- 
size planet around M dwarfs which uses a network of small, identical robotic telescopes 
(optimal bandpass and field of view, telescope aperture, telescope time allocation, observing 
strategy) have been laid out by Nutzman & Charbonneau (2008). That study assumes that 
the well-characterized Mt. Hopkins site in Arizona chosen for the survey guarantees that a) 
the limiting system performance (photometric precision) for transit detection is achieved, and 
b) the requested number of useful observing nights per year is reached. This work presents the 
details of a characterization study carried out at OAVdA, which aims at demonstrating that 
the main site-specific requirements (night-sky brightness, number of useful nights, seeing, 
photometric stability) are satisfied, and OAVdA is therefore an optimal choice for hosting a 
precision photometric search for low-mass, small-size planets around a well-defined sample 
of M dwarfs in the solar neighborhood which utilizes an array of small telescopes. 

The paper is organized as follows. In § 2 we describe the instrumentation presently 
available at OAVdA, which was utilized for the purpose of the feasibility study, and outline 
in § 3 the various steps involved in the dedicated pipeline for the data processing and analysis 
of the photometric data under development. We present in § 4 the observing strategy adopted 
and the main results of the feasibility study in terms of overall site characterization (i.e. the 
achievable photometric precision as a function of atmospheric conditions). We conclude in 
§ 5 by discussing the preparatory steps for a long-term photometric monitoring program 
to characterize the micro-variability characteristics of and search for transiting small-size 
planetary companions to a well-defined sample of low-mass stars, to be carried at the OAVdA 
site in the near future. 
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2. System description 

2.1. Hardware 

For the purpose of the feasibihty study we adopted existing instrumentation available 
at the OAVdA site. Most of the observations were carried out with a 250-mni Maksutov 
telescope (the optical scheme is shown in Fig. [1]), with a focal lenght of 950 mm (f/3.8), on 
a german equatorial mount, and no autoguiding (while not common practice, this does not 
constitute a serious issue for the purpose of our study). The instrument is equipped with 
a front-illuminated Charged Coupled Device (CCD) camera Moravian G2-3200ME (sensor 
KAF 3200ME, area 2184 x 1472 pixeP, pixel area e.SxG.S/im^), and a non-standard As- 
tronomik R filter (see Fig. |2]). As shown in Fig. [3l the sensor reaches a quantum efficiency 
QE ~ 87% at this wavelenght. The CCD camera was operated with a nominal gain of 1 
e~/ADU (in 1 X 1 binning mode), with nominal read noise A^,, = 10 e~ pixel"^ (rms), dark 
current No = 0.1 e~ pixel^^ s~^ at a temperature of — 40°C and a full frame download 
time of ~8.4 s through a USB 2.0 port. The CCD chip is cooled thermoelectrically by a 
two-stages system of Peltier modules which allows to reach a maximum temperature gra- 
dient AT ~ 50 ± 0.1 °C with respect to the ambient temperature. The chip temperature, 
recorded in the fits file headers, is typically stable to better than 0.5 °C. During the period 
of the characterization study (late spring/summer) , the chip was cooled to temperatures of 
— 30°C or — 35°C (while the working temperature is usually set to — 40°C during winter 
time). As shown in Figure HI the typical values of thermal noise for our CCD camera re- 
main well below 1 e~ pixel'^ s^^ as far as the cooling temperature is set to < —25 °C. The 
telescope-CCD configuration has a Field of View (FoV) of 52.10 x 35.11 arcmin^ and a plate 
scale of 1.43"/pixel. 

The OAVdA site also routinely utilizes two Ritchey- Chretien refiector telescopes, one 
400 mm (f/7.64) and one 810 mm (f/7.9), the latter being the primary scientific instrument 
of the observatory. Both the instruments have equatorial open fork mounts. 

The 400 mm telescope will be used full-time together with the 250 mm during the 
next steps following this study, as detailed in § 5. It is coupled to a front-illuminated 
CCD camera Finger Lakes Instrumentation Pro Line PLIOOIE with sensor KAFIOOIE and 
standard BVRI filters. The camera sensor area is 1024 x 1024 pixel^ with pixel area 24 x 24 
/im^. The final configuration is characterized by a FoV of 26.4 x 26.4 arcmin^ with a plate 
scale of 1.55"/pixel. The primary telescope is coupled to a back-illuminated CCD camera 
Finger Lakes Instrumentation Pro Line PL 3041-BB with sensor Fairchild Imaging 3041 and 
standard BVRI filters. The camera sensor area is 2048 x 2048 pixeP, with pixel area 15 x 15 
;Lim^. The system has a FoV of 16.3 x 16.3 arcmin^ with a plate scale of 0.48"/pixel. The 
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detector is used routinely in binning mode 2 x 2, a choice driven by the need to avoid image 
oversamphng, considering the typical seeing conditions at our site (see § l4.4p . While usually 
devoted to other projects, it will be used in the future to perform multi-band photometric 
observations within the context of the necessary follow-up scheme of potentially interesting 
transit candidates. Each observing station is equipped with a pc which is synchronized every 
15 minutes via Internet to the Network Time Protocol (NTP) server of the Italian National 
Institute of Metrological Research. In the header of each frame, time is saved as Julian Date 
( JD) with six significant digits. Table [T] summarizes the specifications of each of the three 
telescope systems. 

The 250 and 400 mm telescopes are located in one of the platforms of the observatory. 
Fig. E] shows a panoramic view of the telescopes location. We evaluated with a theodolite the 
minimum sky altitude above the horizon, as a function of the azimuth, that can be reached 
with a telescope placed in the middle of the platform before encountering any obstacle 
limiting the observations (mountains, platform walls, dome). The altitude-azimuth profile 
is shown in Fig. [6l The most prominent obstacle is represented by the dome (with a 7-m 
diameter) hosting the primary telescope. The dome (not shown in Fig. \5^ is located next 
to the platform. It covers ~ 20 — 24 deg in elevation, in the azimuth range ~290-320 deg. 
This does not represent a serious constraint, however, as per standard procedure we observe 
targets until they reach a minimum altitude of 30 deg above the horizon. 



2.2. Telescope control system 

The 400 and 810 mm telescopes are managed by the control and data acquisition system 
CompactRIO which is based on National Instruments (NI) technology. The system is pow- 
ered by a software with a graphical interface developed with the NI Lab VIEW programming 
tools and it allows the remote control of the telescopes (the pointing is both automatic and 
in manual jog). The algorithms integrated in the control system were entirely developed by 
members of the OAVdA reasearch staff. The system is able to dialog, through the use of the 
LX200 protocol, with the most popular softwares used to control astronomical instruments. 



2.3. Seeing monitoring 

During the period of the characterization study, each night of observation we monitored 
and recorded the seeing using the Differential Image Motion Monitor (DIMM) technique 
(Sarazin & Roddier 1990), which minimizes the effects induced by the instrumentation and 



provides an accurate estimate of the level of the atmospheric turbulence. This technique 
requires the use of a two-hole mask (Shack-Hartmann mask) in front of the telescope aperture 
so that, slightly defocusing, two distinct and close profiles of one bright star are formed. Fig. [7] 
shows a freeze-frame taken during a typical seeing monitoring session illustrating the DIMM 
procedure. A seeing estimate is then obtained measuring the relative motion of the two 
stellar profiles over time scales of ms. To monitor the seeing we used a second 250 mm f/3.8 
Maksutov telescope (shown in Fig. [5]) equipped with a Shack-Hartmann mask and a SBIG 
STV digital integrating video camera which is able to apply in real-time the DIMM algorithm 
with high time sampling (typically, 5 measurements per second). Data are automatically 
saved on a PC as ASCII files. No filter was used to perform these measurements (possibly 
resulting in an overestimate of the seeing). 



3. Data processing and analysis 

A fundamental element of the characterization study has been the development and 
testing of the automatic pipeline TEEPEE (Transiting ExoplanEts PipElinE) for the re- 
duction and processing of the photometric data. TEEPEE is a software package written 
in IDLO, which utilizes publicly available software from the Astronomy Users' Library U as 
well as external contributed FORTRAN routines. Two routines, aper and getpsf , are IDL 
versions of original IRAF routines |j. The original portion of software is the one devoted to 
automatically perform ensemble differential aperture photometry on an user-specified stellar 
target. The various tasks performed by TEEPEE and the different routines used are sum- 
marized in the diagram shown in Fig. [HI TEEPEE is organized in three main, sequential 
blocks, whose functionalities are described below. 



3.1. Image calibration 

In the package imgcalib, every raw FITS image acquired by the CCD (each 6.1 Mb 
in size) is initially dark-subtracted and fiat-fielded. A master-dark image is created every 



^IDL is a commercial programming language and environment by ITT Visual Information Solu- 
tions, http : //www . ittvis . com/idl/ 



http : //idlastro . gsf c . nasa.gov/contents . html 



^IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Asso- 
ciation of Universities for Research in Astronomy (AURA) under cooperative agreement with the National 
Science Foundation. 
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month from the median values of each pixel of 70 dark images having the same exposure time 
of the science imageq^. The flats are prepared from the nightly averages of ~15 high-S/N 
exposures at twilight, diffused by a white opal of polistyrene. 



3.2. Photometric & astrometric processing 

The first task of the pipeline package photomak is the search for the photocenters in 
each image. This is performed automatically by the subroutine detectcenter, which selects 
each pixel with an ADU value greater than 3 times the rms of the sky background. The list 
of X and y coordinates of the photocenters are used as input for the subroutine psf , which 
calculates the average FWHM of the point-spread function (PSF) of each detected star in 
every image. Then, the aper subroutine is applied to perform single-aperture photometry 
on each detected object, with the radius of the user-provided circular aperture typically set 
to 2 — 3 times the average value of the FWHM. This routine counts the ADU values for 
each pixel in the circular aperture, while partial pixels are handled by numerical integration. 
Once the photocenter coordinates and the instrumental magnitudes are listed, the daomatch 
subroutine alignes each image to a reference frame (usually, the first image of each series) 
by performing a triangulation. If in one image the stars found are less than 20% of those 
detected in the reference frame or the alignment procedure fails (mostly because of high 
variability in the sky transparency), then the image is automatically discarded. 



3.3. Data analysis 

The heart of TEEPEE, for the purpose of this study, is the package analysis. It 
performs those operations which are necessary to eliminate, when possible, all the systematics 
which cause the degradation of the photometric precision, and outputs the photometric light 
curves for stars detected in a every image of a field. At first, by using the subroutine 
starlist only the stars that have been detected in every frame are selected and added 
to the final list. All objects that, due to telescope drifts exited the CCD FoV, or are too 
weak (i.e. below 3cr of the sky background level) in one or more frames are discarded. The 
typical drift of the 250 mm telescope is nearly 10 pixel hr~^ and 20 pixel hr~^ in the X 



^'^The number of dark images ensures minimization of the noise for the purpose of achieving a goal of a few 
milh-mag photometric precision. We are aware of the fact that taking dark images during each observing 
session might improve to some degree the quahty of the cahbration. The protocol will be changed in the 
future, by taking for example dark frames both at the beginning and at the end of the observing session 
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and Y direction, respectively. Considering, for example, the fields of WASP-3 and HAT-P-7 
which we monitored for this study (see § 14.31) . on average ~ 15% of the stars appearing 
in the first frame of a series are lost because of drifts or veils, resulting in some 550-600 
good stars detected on average in every frame during a night. Finally, the corrections due 
to atmospheric extinction are applied. The equation giving the measured magnitudes as a 
function of the magnitude above the atmosphere is: 

mj = mo + kx AM, (1) 

where k denotes the atmospheric extinction coefficient, rrij the instrumental magnitude, 
rrio the magnitude above the atmosphere and AM the air-mass {AM = 1 when the object is 
located at the local zenith). The standard procedure requires that slope k and intercept rrio 
be evaluated every time from observations of photometric standard stars, such as Landolt 
stars. Given the purely differential approach to the photometric measurements in this study, 
we proceeded to calculate the extinction coefficient using all the n reference stars chosen in 
the field of each target. This is done by solving a global system of nxNi linear equations 
of the form of Eq. [T] through the IDL routines SVDC and SVSOL, where iVj is the number of 
good frames of the series. The solution is represented by a global value for k and Ni values 
for mo- As comparison stars we selected the n brightest stars detected in every frame with 
mj < 13. For the purpose of this study, the i?-band measurements of k are used to gauge 
the level of sky transparency (see § 14.11) . No dedicated device (e.g., CAVEX, Sanchez et al. 
2007) for the monitoring of extinction was used. 

The rms of each new light curve obtained after the air-mass correction provides quick 
information about the stability of the atmospheric transparency during the night, but more 
accurate light curves are provided by the ensemble differential photometry. For each frame 
i, we used as the reference magnitude m^^r the average magnitude of the n reference stars 



</ = ^^^^^^. (2) 



iRrif is then subtracted to the magnitude of the user-defined target Mtarget, obtaining the 
difference Am* = m^^^ — m\^^g^^. The procedure is iteratively repeated for all the reference 
stars too, using as new references the remaining n — 1 stars. We thus obtain a series of 
differential light-curves for all the stars selected, which can be inspected for variability. The 
few bright variables identified with this approach (see § 14. 6p are not discarded at this stage, 
as their impact on the photometric precision is minimized by the selection of dozens of 
comparison stars. 
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As discussed, the version 1.0 of the TEEPEE pipehne can perform automatically several 
basic tasks. These have been tested with promising results, as described in the following 
sections. However, the pipeline is still far from being completed. Further refinements will 
guarantee improved functionality and a higher degree of robustness. In particular, a pho- 
tometric survey finalized to the detection of planetary transits must tak e into account th e 



effects of 'red noise', which can strongly affect the sensitivity thresholds (jPont et al.ll2006l ). 
This feature will be included in future upgrades of the software package in preparation of 
the official start of the photometric survey in the near future. The pipeline will grow further 
to implement algorithms for the reduction of the noise and systematic variations in the light 
curves (detrending), such as the Trend Filtering Algorithm (TFA, Kovacs et al. 2005) or 
SysRem (Tamuz et al. 2005), as well as algorithms for the detection of pe riodic transit- 



like s ignals, such as the widely used Box-fitting Least Squares (BLS) method (JKovacs et al. 



2002|). 



4. Site characterization study: results 

The main OAVdA site characterization campaign was carried out during the period 
May- August 2009. As in any search for adequate sites for new observatories, we focused 
on the determination of observing conditions such as night-sky brightness, number of useful 
nights, photometric precision, and seeing. 



4.1. Night-sky brightness and transparency 

We have utilized measurements of the zenithal night-sky brightness collected between 
May 2007 and July 2007 during seven moonless nights in collaboration with the Italian 
Agency for Environmental Protection (ARPA), as part of an extended monitoring program 
of the light pollution in the Aosta Valley. The data were obtained using the 400 mm re- 
flector telescope and a CCD camera equipped with a l^-band filter. All measurements were 
calibrated using standard reference stars. The y-band night-sky brightness ranged between 
21.12 and 21.61 mag arcsec"^, with a mean and standard deviation of 21.29 ± 0.16 mag 
arcsec^^. Th ese numbers are ty pical of those of very good, very dark observing sites (e.g.. 



see Table 4 in lMoles et al.l (|2009[ )). During the same period, the \^-band extinction coefficient 
ranged from 0.116 mag to 0.420 mag, with a mean and standard deviation of 0.26±0.13 mag. 
This value compares reasonably well with the l^-band k estimates during summer available 
in the literature for other sites (e.g., Sanchez et al. 2007, and references therein). The mean 
and standard deviation of our own i?-band k value derived within the context of the site 
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characterization campaign are 0.165 ± 0.090 mag. They are in good agreement with the 
analytical estimate, 0.146 mag, from the extinction curve formula in Sanchez et al. (2007) 
based on the \^-band value of k reported above. The typical conditions of sky transparency 
are thus in line with what to expect from a good site, although the spread around the mean 
conditions is relatively large. Given the limited time span of the observations, we cannot 
comment on seasonal variations of the extinction at the OAVdA site, nor are we in a position 
to investigate in detail the relative role of various contributors (e.g., dust and other aerosol 
particles, ozone). 



4.2. Fraction of useful nights 

In order to estimate the number of useful nights/yr, we relied upon on a 3-yr baseline 
of data (from 2006 to 2008) collected by a weather station located next to the Observatory, 
while for the year 2009 we also used our own internal logs. Over this time span, fully 
clear nights conditions are statistically realized 142 night/yr (39%), while an additional 67 
nights/year (18%) can still be considered useful, assuming a typical cloud cover of < 50% of 
the useful observing time (conservative upper limit). As for the year 2009, during which we 
carried out the feasibility study, we registered similar conditions: uninterrupted observations 
could be carried out 39.5% of time, while for an additional 22% of the nights the cloud cover 
did not exceed 50% of the useful observing time. These figures compare reasonably well with 
other observatory sites in Europe. 



4.3. Photometric performance 

In order to estimate the achievable photometric precision at the OAVdA site, given 
the above mentioned instrumental setup, the centerpiece of the site characterization study 
consisted in the monitor i ng of three stella r fields centered on the transiting pl anet hosts 



Gliese 436 (IButler et al.ll2004[ ). WASP-3 flPoUacco et al.M2007h . and HAT-P-7 (IPal et al. 



20081 ). The targets were observed during both in- and out-of-transit phases during the 
period May-August 2009. All observations were performed with the CCD set up in 1 x 1 
binning mode and at the focus of the telescope, and the exposure times were chosen at 
the beginning of each session (without being subsequently modified) in order to guarantee 
an optimum signal-to-noise ratio {S/N) for the target, while avoiding saturation. Typical 
exposure times were 22 s for Gliese 436 (y=10.68 mag), the brightest target at R band, 
and ranged between 35 and 60 s for both WASP-3 (1^=10.48 mag) and HAT-P-7 (^=10.46 
mag), given variable conditions of air mass and sky transparency. Fig. IH] shows an example 



-13- 



of the quality of the images we collected. The whole database for this characterization study 
comprises ~10,300 good images, corresponding to ~ 63 GB of data. 

We show in Fig. (TU] the scatter in our differential photometry as a function of apparent 
magnitude in the field of the star WASP-3, for a night representative o f the typical perfor- 



manc e achieved. The results are compared to a theoretical noise estimate (JNutzman fc Charbonneau 



20081 ) which includes contributions from the object (Poisson noise), sky, readout noise, and 



scintillation: 



(^ph = j;^ (3) 

In Eq. [3], A'^ is the number of received photons from the star during the exposure, Upix 
is the number of pixels in the photometric aperture, A^^ is the average number of photons 
per pixel of the sky background (we used a number typical of the average sky conditions 
registered over the duration of our campaign). No is the dark current contribution in e~ 
pixel"^ (from § 12.11) . Nr is the rms of the CCD read noise in e~ pixel"^ (also from § 12. ip . 
and (Jscint is the noise term for the atmospheric scintillation. The expression used for agdnt 
is (Dravins et al. 1998): 

'-' 0.09^^expr-^V (4) 



where AM is the air-mass (we adopted a mean value of 1.5), D is the telescope aperture 
in cm, t is the exposure time in s, and h is the height in km of the observing site above the 
sea level. 

For bright stars, the dominant sources of error are the scintillation and shot (Poisson) 
noise from the target star itself. Notice the excellent agreement with the expected scatter 
values for magnitudes i? < 13 mag. Among the most likely sources of additional scatter 
in the observations with respect to the theoretical predictions we list guiding errors and 
possible contributions from centering errors during aperture photometry. For the brightest 
objects we achieve a nightly rms precision of ~ 0.003 mag (during some nights the recorded 
scatter for these stars was even lower). For fainter stars, sky photon noise is the dominant 
source of scatter. Some outliers can be seen in the dat aset. One of them c orresponds to a 



newly discovered variable (probable 6 Scuti-type) star (JDamasso et al.ll2010l ). 



To evaluate the photometric precision of each night of observation we adopted the 
mean rms calculated for stars with magnitudes i? < 13 mag. In Fig. [11] we show how the 
nightly-averaged rms varied during the period of observation. Different symbols are used to 
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distinguish the various star fields that were monitored, as described in the figure caption. 
The median rms over the whole duration of the characterization study is 0.006 mag for 
magnitudes brighter than 13 mag. The results shown in Fig. (TU] can then be considered 
representative of a typical night at our site. Equally important is the fact that the night-to- 
night variations are ~ 0.002 mag (disregarding the few nights with large, > 0.01 mag scatter 
due to poor weather conditions). These results demonstrate that it is possible to achieve 
competitive precision at the OAVdA site over an extended period of time, a particularly 
encouraging result given the still non-optimized hardware system used during this study. 



4.4. Seeing measurements 

As described in § 12.31 seeing monitoring has been carried out in conjunction with the 
photometric measurements over the whole duration period of the site characterization study. 
For the purpose of this work, its role has been quantified in terms of its median value over 
the timespan of the measurements, its typical stability during a night, and its correlation 
with the photometric performance achieved. We note that, given the limited time duration 
of the measurements, we cannot comment in a statistically significant fashion on possible 
seasonal seeing variations, or in detail on its possible dependence on other meteorological 
parameters, such as wind speed and direction. 

Fig. [12] shows the distribution of seeing values as a function of the night of observation. 
The median seeing turns out to be ~ 1.7" over the whole period. In most nights, the seeing 
appeared rather stable, typically remaining within 20% or so of the reference median value. 

Indeed, the overall result is that the OAVdA site does not appear competitive with 
respect to other observing sites in the world (see, e.g.. Table 5 in Moles et al. 2009). However, 
as illustrated in Fig. [13] the photometric precision is essentially seeing independent, for values 
up to 3". On the basis of this result we can conclude that the typical seeing measured at the 
OAVdA site does not represent a limiting factor for the precision of the stellar differential 
photometry (with the instrumental set up we utilized), for magnitudes i? < 13 mag. This is 
in accordance with the promising results found applying the telescope defocusing technique 
for high-precision photometry of transiting extrasolar planets (i.e. Southworth et al., 2009a, 
2009b, 2009c), which can be carried out for the brightest stars even with small-size telescopes. 
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4.5. Sample light curves of known transiting systems 

As examples of the quality of the photometric data achieved with the present instru- 
mental setup, we present two light curves showing the transits of the extrasolar planets 
WASP-3b (Fig. [14]) and HAT-P-7b (Fig. [15]). As for Gliese 436b (whose data are not shown 
here), unfavourable wheater conditions prevented us from observing a full transit during the 
timespan of the characterization study. 

WASP-3 (M^ = 1.24M0, R^ = 1.31 i?©), which hosts an Mp = 1.76 Mj, Rp = 1.31 Rj 
planet with a period P = 1.846834 days, was monitored for 19 nights. The rather good phase 
coverage allowed us to observe four transits of WASP-3b. Fig. [TH shows the full transit 
observed during the night of UT 28 July 2009. Superposed to the normalized flux is the 
best-flt curve, obtained as solution of a non-linear least square probler n using the Levenaberg- 



Marquardt algorithm and a merit function based on the formalism of iMandel fc Agoll (l2002[ ) 
with adjustable quantities the inclination angle i, the ratio of radii r = Rp/R^^, and the 
impact parameter b = {a/R^,)cosi (where a/R^, is the orbital semi-major axis expressed 
in units of the stellar radius). We assumed a quadratic limb darkening law with i?-band 
coeflicients mi=0.2357 and ^2=0.3779 (Claret 2000), appropriate for the stellar parameters 
in PoUacco et al. (2008). The sub-panel of Fig. UM shows the post-fit residuals, that exhibit 
an rms of 0.002 mag. The relatively large transit depth (0.013 mag) allowed us to derive 
formally precise values for the system parameters. The best-fit values and their associated 
1 — cr uncertainties derived from the covariance matrix of the solution are as follows: b = 
0.665 ± 0.003, i = 81.24 ± 0.06 deg, and r = 0.1091 ± 0.0006. There appears to be some 
discrepancy (a less central transit) with respect to the values reported in the discovery paper 
of PoUacco et al. (2008) and in the work of Gibson et al. (2008). However, the results of this 
exercise are shown here only for illustrative purposes. Our results are undoubtedly affected 
by the non-optimal choice of the limb-darkening parameters (we recall observations were 
taken with a non-standard R filter), by the lack of a proper treatment of red noise, and by 
the almost complete absence of data gathered before ingress. 

The star HAT-P-7 (M^ = 1.47 Mq, R^ = 1.84 Rq) was selected as a challenging target 
for our test because its transiting planet {Mp = 1.8 Mj, Rp = 1.42 Rj, P = 2.2047298 days) 
produces a photometric transit characterized by a shallow depth (~0.007 mag), which is 
close to the detection limit for our present instrumental and data reduction and analysis 
setup. HAT-P-7 is located in the field of view of the Kepler space telescope, that monitored 
the star for t en days during its c ommissioning phase, deriving an exquisite-quality optical 



phase curve (JBorucki et al.ll2009l ). In Fig. [15] we show the transit observed at OAVdA the 
night of UT 18 July 2009. We could not observe a complete transit (the planet's egress is 
barely visible due to the approaching sunrise which forced us to interrupt the exposures). 
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and during the night the sky transparency varied significant due to the presence of veils. 
However, the transit was clearly detected, and by using an aggressive light-curve analysis, 
including iterative a clipping for the removal of outliers, we managed to 'clean' the dataset 
so as to perform an analogous exercise to that carried out in the case of WASP-3. The 
light-curve fit, assuming a quadratic limb darkening law with i?-band coefficients mi=0.3980 
and M2=0.2071, appropriate for the stellar parameters in Pal et al. (2008), resulted in the 
following solution: h = 0.321±0.006, i = 85.93±0.08 deg, and r = 0.064±0.002. In this case, 
the main difference is that of a value of r which is about 15% smaller than that reported in the 
literature. This is due to the shallower transit depth (~ 5 mmag) recovered after we applied 
the procedure to remove outliers. Again, we stress the illustrative purpose of the orbital fit 
carried out, the main aim here being simply that of showing that the transit could clearly be 
detected. This represents a promising result in the perspective of our long-term photometric 
monitoring program to be started in the near future with upgraded instrumentation and an 
improved data reduction and analysis pipeline. 



4.6. New variable stars 

A natural byproduct of high-precision, high-cadence photometry collected for hundreds 
of stars over a timespan of several weeks, as is the case for the WASP-3 and HAT-P-7 
stellar fields monitored in the context of our site characterization study, is the discovery 
of variable sources. A handful of previously unknown variable stars with relatively large 
amplitude variations (AV" > 0.05) were uncovered during this study. Two of them are 
presented and discussed in Damasso et al. (2010). Here we report the detection of four 
new binary/multiple-star systems showing eclipses. The fields centered on the variables are 
shown in Fig. [121 while the relevant data are summarized in Table [21 

We searched the literature and stellar catalogues in order to find prio r variability 



i nforrn atio n on the four objects. In particular, we looked at the works of [Street et al. 



( |2007[ ) and [Norton et al.l (120071 ). and we searched the General Catalogue of Variable Stars 



(GCVS; http://www.sai.msu.su/groups/cluster/gcvs/gcvs/), the New Catalogue of 



Suspected Variable Stars (http://www.sai.msu.su/groups/cluster/gcvs/gcvs/nsv/ ), the 



VSX and ASAS databases (http : //www . aavso . org/vsx/'and'http : //www . astrouw . edu . pl/asas/?page= 
and the VizieR database (http://vizier.u-strasbg.fr/). We failed to identify counter- 
parts of our objects, so that we can report them here as new discoveries. 

The first two targets in Table M (GSC2.3 N208000215 and GSC2.3 N2JH035417) are 
identified as eclipsing binary systems, with orbital periods P = 2.02065 ± 0.00004 days and 
P = 0.505167 ±0.000002 days, respectively, based on a Lomb-Scargle periodogram analysis. 
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In the two panels of Fig. [T7]we show the phase- folded light curves for these two systems. For 
the first binary system (upper panel of Fig. [T71) we derived the ephemeris of the secondary 
minimum, that was observed at the predicted epochs even if its depth is shallow (~0.005 
mag) and very close to the detection limit for our instrumentation. 

For the last two targets in Table [21 we observed in their light curves only one complete 
minimum (GSC2.3 N2JH000428) and a partial minimum (GSC2.3 N2JH066192), as showed 
in the upper and lower panels of Fig. [TS] respectively. Because we observed these minima just 
one time, without detecting any other evident variability in other data, we can not estimate 
the orbital periods of the systems and we can only classify them generally as eclipsing binary 
stars. 



5. Concluding remarks 

We have carried out a site characterization campaign at the Astronomical Observatory 
of the Autonomous Region of the Aosta Valley (OAVdA), in order to demonstrate that 
OAVdA can be the location for a long-term photometric survey with small telescopes aimed 
at detecting transiting rocky planets around nearby cool dwarf stars. 

We focused the attention on those site-dependent factors that can have the largest 
impact on the ultimately achievable precision of the photometric measurements, such as 
seeing, extinction and night-sky brightness. We then correlated them with the quality of the 
photometric measurements of selected target fields, monitored with 20-40 cm class telescopes 
available on site, and analyzed using standard techniques of differential aperture photometry 
using an ad hoc developed data processing and analysis pipeline. Relevant data were collected 
in situ over a period of ~ 4 months, and complemented by additional meteorological and 
photometric datasets covering a timespan of > 3 years. The main findings of this study 
(partly presented already in Damasso et al. 2009) can be summarized as follows: 

• The measured zenithal \^-band night-sky brightness is typical of that of very good, 
very dark observing sites. The extinction registered at V band is not dissimilar from 
that of other sites (for the same season) . 

• The median seeing over the period of in situ observations is found to be ~ 1.7". Given 
the limited duration of the observations, we did not probe any possible seeing seasonal 
patterns, or the details of its possible dependence on other meteorological parameters, 
such as wind speed and direction. Moreover, our data show that the seeing at the 
observatory was reasonably stable during most of the nights. 
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• The fraction of fully clear nights per year amounts to 39%, while the total of useful 
nights increase to 57% assuming a typical cloud cover of not more than 50% of the 
night. There is not yet enough statistics to determine a clear value for the fraction of 
photometric nights per year. 

• During 24 good observing nights over a period of over 3 months the median photo- 
metric precision was 0.006 mag for stars with magnitude i? < 13 mag. The typical 
nightly photometric precision appears to be uncorrelated with the seeing, whose typical 
value (in principle not competitive with other highly reputed observing sites around 
the world) turns out not to be a limiting factor for achievement of the photometric 
performance required for detection of planetary transits. 

• The observation of known transiting planets (those in particular of WASP-3b and HAT- 
P-7b) allowed us to show that our differential aperture photometry approach applied 
to data collected with a small, 250-mm telescope, achieves best-case rms precision of 
0.002 — 0.003 mag for the brightest stars. 

• A byproduct of the significant amount of photometric data collected in the stellar fields 
of WASP-3 and HAT-P-7 was the identification of a handful of new variable stars, four 
of which are presented and discussed here. 

The results of the site testing campaign show conclusively that OAVdA is a suitable 
choice for hosting a long-term photometric survey for transiting planets around cool stars 
in the solar neighborhood. A preliminary step towards the definition of such a long-term 
photometric transit search program is the 'pilot study' which we have begun earlier this year. 
We are devoting one of the 250 mm telescopes and the 400 mm instrument to the J-band 
photometric monitoring of ~ 50 M dwarfs, within ~ 25 pc from the Sun, with high-precision 



paral lax estimates determined within the context of the TOPP program (jSmart et al.l 12007 



20101 ). This study, with a targeted duration of less than 1 year, aims at obtaining uniform 
photometric monitoring of all targets (on the order of at least a few weeks of data per target), 
with the goal of providing benchmarks for the photometric micro- variability characteristics 
of M dwarfs as a function of their spectral sub-type (the stars in our sample span the range 
MO through M9), in relation to other, primarily spectroscopic, activity indexes (e.g., based 
on Hq,, Ca II lines emission). While the pilot study will likely not detect any planetary 
transits (given the limited number of objects in the target sample), the large photometric 
dataset collected will help us placing useful reference points for improved understanding of 
the impact on transiting planets detectability of stellar activity- related phenomena (e.g., 
spot-induced rotational modulation), and the prospects for taking them into account via 
detailed modeling procedures. 
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In parallel to the science-driven activities of the pilot study, we are making progress 
in the re-organization of the scientific terrace for the purpose of hosting multiple identical 
400- mm robotic telescope systems (at least five), the start of the operations of the first pair 
being envisaged in connection with the end of the pilot study (Fall 2010). The hardware and 
software requirements for the control system for data taking and storage have been clearly 
defined. Finally, as mentioned in § 13.31 necessary upgrades to the TEEPEE automated data 
reduction and analysis pipeline are being considered, in order to increase its functionality 
and robustness, ahead of the official start of our Western Italian Alps photometric transit 
search for planets around M dwarfs. 
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Table 1: Summary of the main characteristics of the telescope systems described in § 2. 



Telescope 



CCD camera 



Resulting configuration 



Optical scheme 



Aperture Focal ratio Sensor area Pixel area FoV 

(mm) (pixel) (/^^i) (arcmin) 



Plate scale 
("/pixel) 



to 
to 



Reflector Maksutov 
Reflector Ritchey- Chretien 
Reflector Ritchey-Chretien 



250 


f/3.80 


400 


f/7.64 


810 


f/7.90 



2184 X 1472 6.8 x 6.8 52.10 x 35.11 1.43 (binning 1x1) 
1024 X 1024 24 x 24 26.4 x 26.4 1.55 (binning 1x1) 
2048 X 2048 15 x 15 16.3 x 16.3 0.48 (binning 1x1) 



Table 2: Variable stars discovered in the fields of the extrasolar planet hosting stars WASP-3 and HAT-P-7. R 
magnitudes are derived from r' magnitudes of the CMC14 catalogue using the relation R = r' — 0.22 derived by 
Dvmock fc Milesl fcoooh 



Object ID 



Eq. Coord. 



Magnitude 



(GSC2.3 catalogue) 



RA, DEC (J2000) 



B 


V 


R 


J 


1^ 00 


12.34±0.19 


12.09±0.18 


11.622 


10.962±0.02 


10.662±0.023 


- 


- 


12.797 


12.365±0.022 


12.225±0.024 


13.34±0.30 


14.10±0.44 


12.95 


11.998±0.020 


11.053±0.019 


10.11±0.04 


8.92±0.02 


- 


6.857±0.018 


6.205±0.016 



N208000215 
N2JH035417 
N2JH000428 
N2JH066192 



18:35:23.210 +35:14:45.43 
19:27:49.596 +48:20:47.90 
19:29:55.081 +47:50:04.41 
19:30:08.505 +47:49:31.16 
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Fig. 1. — Optical scheme of the Maksutov reflector telescope (aperture 250 mm) used for 
the characterization study (from|http://lnx.costruzioniottichezen.com|). 
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Fig. 2. — Transmission curve for the set of filters coupled to the CCD in the 250nini telescope. 
For the site testing we used the R filter. The corresponding transmission curve is showed in 
red. 
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Fig. 3. — Nominal response curve of the sensor KAF-3200ME integrated in the CCD camera 
used in this study (solid line). 
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Fig. 4. — Thermal noise as a function of temperature for the sensor KAF-3200ME (pixel 
size 6.8 X 6.8 /im^) integrated in the CCD camera used in this study (green solid line). The 
performances of other off-the-shelf devices are shown for comparison. 



-28- 




Fig. 5. — The OAVdA scientific platform from wliicli tlie site cliaracterization study ob- 
servations were carried out, and wliere an array of telescopes will be installed for a future 
long-term photometric survey. Numbers 1 and 3 identify the two 250 mm Maksutov tele- 
scopes used, respectively, to monitor the seeing with the DIMM technique and to acquire 
photometric data of fields with stars hosting transiting extrasolar planets. In the middle 
(number 2) is the 400 mm reflector telescope. 
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Fig. 6. — Minimum sky altitude useful for observations as a function of the azimuth, before 
incurring into obstacles (observatory dome, platform walls, and mountains). The profile was 
determined with a theodolite placed in the center of the OAVdA platform which hosts the 
telescopes used in this study. 
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Fig. 7. — Freeze-frame from the SBIG-STV digital integrating video camera coupled to a 250 
mm Maksutov telescope to monitor the seeing with the DIMM method. In the upper left 
corner the two images of the star used to estimate the seeing are represented, as produced 
placing a two-hole Shack-Hartmann mask in front of the telescope aperture. 
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Fig. 8. — Flowchart summarizing the main tasks (and corresponding routines) performed by 
the version 1.0 of TEEPEE, the pipehne for data reduction and processing developed during 
the study presented in this paper. More details are given in the text. 
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Fig. 9. — A typical FITS image obtained with the 250 mm telescope. In this example the 
field of the star HAT-P-7 is shown, indicated with a cross (Exposure time: 60 s). 
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Fig. 10. — Photometric errors (rms) vs. instrumental mean magnitude in the field of the 
planet-hosting star WASP-3 during a good observing night (and for an exposure time of 35 
s). The various contributions to the expected photometric noise, according to Eq. |3l are also 
shown. 
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Fig. 11. — The nightly-averaged photometric precision, calculated for stars with magnitudes 
< 13 mag, is plotted for each night of observations. Different symbols are used to distinguish 
the three stellar fields monitored during the characterization study: a plus denotes the field 
of the star Gliese 436, an asterisk is used for WASP-3, while squares correspond to HAT-P-7. 
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Fig. 12. — The distribution of seeing measurements over the whole period of observation. 
The small black stars show 7864 individual data points (many of them overlapping due to 
the seeing values having been approximated to the first significant digit), each corresponding 
to an average seeing value over a 1 min interval, rounded to 0.1". The large filled circles 
indicate the median seeing value for each night. 
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Fig. 13. — The comparison of the nightly- averaged photometric precision to the corre- 
sponding mean seeing does not reveal an evident correlation between these two variables 
for FWHM values < 3", given our intrumental set up. Symbols are as in Fig. [11] 
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Fig. 14.— Top: transit light-curve of WASP-3b during the night of UT 28 July 2009. The 
time of transit center occurred at T^ = 2455041.411 HJD. The 358 35-sec exposures were 
collected over 4.6 hrs. The best-fit curve is superposed. Bottom: the post-fit residuals 
exhibit an rms of 0.0021 mag 
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Fig. 15. — Normalized light-curve of HAT-P-7 showing the shallow transit of its Jupiter-mass 
companion HAT-P-7b (UT 18 July 2009, 260 60-sec exposures collected in 4.6 hrs). The 
time of transit center occurred at Tc = 2455031.5222 BJD. The best-fit curve is superposed. 
Bottom: the post-fit residuals exhibit an rms of 0.0018 mag. 
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Fig. 16. — The new variable stars discovered at OAVdA and presented in this work. A: 
GSC2.3 N208000215. B: GSC2.3 N2JH035417. C: GSC2.3 N2JH000428. D: GSC2.3 
N2JH066192. 
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Fig. 17. — Phase-folded light curves of two eclipsing binary systems discovered during the 
feasibility study. Upper panel: the estimated orbital period is P = 2.02065 ± 0.00004 d and 
phase=0.0 corresponds to HJD=2455041.50975 + P x E. Lower panel: the estimated orbital 
period is P = 0.505167 ± 0.000002 d and phase=0 corresponds to HJD=2455031. 381135 + 
P X E. In the upper and lower panels, there are 1510 and 1443 data points, respectively, 
taken over 11 and 6 days, respectively. 
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Fig. 18. — Light curves of the third and fourth object in Table |2] (field of the star HAT-P-7) 
showing the occurrence of eclipses likely due (based on the observed depths) to a stellar 
companion. Both signals are seen only once in our data. The minimum in the upper panel 
was detected during an observing session aimed at gathering more data for the eclipsing 
binary system showed in the lower panel of Fig. [T71 In the upper and lower panels, there are 
317 and 313 data points, respectively, taken over a single night. 



